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ABSTRACT 

r*j \ We present a 1200- /im image of the Great Observatories Origin Deep Survey North (GOODS- 

Q^. N) field, obtained with the Max Planck Millimeter Bolometer array (MAMBO) on the IRAM 

30-m telescope. The survey covers a contiguous area of 287 arcmin 2 to a near-uniform noise 
level of ^0.7 mJy beam -1 . After Bayesian flux deboosting, a total of 30 sources are recovered 
(^3.5-c). An optimal combination of our 1200-//m data and an existing 850-/im image from 
' the Submillimetre Common-User Bolometer Array (SCUBA) yielded 33 sources (^4-cr). We 

combine our GOODS-N sample with those obtained in the Lockman Hole and ELAIS N2 
fields (Scott et al. 2002; Greve et al. 2004) in order to explore the degree of overlap between 
1200-/im- and 850-/im-selected galaxies (hereafter SMGs), finding no significant difference 
between their 5850 i um/'S , i200 i uin distributions. However, a noise-weighted stacking analysis 
yields a significant detection of the 1200-/im-blank SCUBA sources, 5850 i uin/<S , i200/«n = 
\ 3.8 ± 0.4, whereas no significant 850-/im signal is found for the 850-/im-blank MAMBO 

sources (S850 jU m/' S 'i20u>m = 0.7±0.3). The hypothesis that the Sana^m/ Simoom distribution 
of SCUBA sources is also representative of the MAMBO population is rejected at the ~4- 
a level, via Monte Carlo simulations. Therefore, although the populations overlap, galaxies 
selected at 850 and 1200/im are different, and there is compelling evidence for a significant 
| 1200-/im-detected population which is not recovered at 850 /im. These are submm drop-outs 

(SDOs), with 5 f 85o^m/S'i2ooMm = 0.7 — 1.7, requiring very cold dust or unusual spectral 
energy distributions (Td ~ 10 K; (3 ~ 1), unless SDOs reside beyond the redshift range 
observed for radio-identified SMGs, i.e. at z > 4. 
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1 INTRODUCTION 

A decade ago, surveys undertaken at submillimetre (submm; 850- 
/im) wavelengths with SCUBA (Holland et al. 1999) transformed 
the accepted understanding of galaxy formation and evolution by 
revealing an unexpected yet significant population of dusty, far- 
infrared-luminous starburst galaxies (Smail, Ivison & Blain 1997; 
Hughes et al. 1998; Barger et al. 1998). Since then, great strides 
have been made in understanding the nature of submm- selected 
galaxies - see Blain et al. 2002 for a review - largely facilitated by 
accurately pinpointing SMGs in deep, high-resolution radio maps 
(e.g. Ivison et al. 1998a, 2000, 2002; Smail et al. 2000) which helps 
overcome the coarse resolution of SCUBA (FWHM ~ 15 arcsec). 
The majority (60-80 per cent) of bright (Sssoum ^ 5mJy) SMGs 
have been identified in this way and one of the most important 



steps forward was the acquisition of about 100 optical/near-infrared 
spectra of such radio-identified SMGs. This placed them typically 
in the redshift range z ~ 1—3, with a median of 2.3 (Chapman et al. 
2003, 2005). However, while submm surveys are equally sensitive 
to sources of a fixed luminosity in the redshift range z ~ 1 — 8, due 
to the negative fc-correction offsetting the cosmic dimming (Blain 
& Longair 1993), the detection rate of SMGs in even the deepest 
radio surveys drops off rapidly beyond z ~ 3.5. Thus, the spectro- 
scopic requirement for a robust radio detection biases the known 
population to z g 3.5. This led to the suggestion that the 20^10 
per cent of the bright, radio-blank SMGs (1.4 GHz flux density, 
5i.4GHz S 20/iIy) could be the high-redshift tail of the popula- 
tion (z > 4 - see Carilli & Yun 1999; Ivison et al. 2002; Eales et 
al. 2003; Aretxaga et al. 2003, 2007), although using photometric 
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redshifts Pope et al. (2006) argued that < 14 per cent of the bright 
SMG population is at z > 4. 

Perhaps the strongest candidates for z > 4 SMGs have come 
from submm/mm interferometry, with the IRAM Plateau de Bure 
Interferometer (PdBI) and the Submillimeter Array (SMA), of a 
few very bright SMGs which - by virtue of their radio, submm/mm 
and near-infrared properties - were deemed to reside at very high 
redshift (Dannerbauer et al. 2002; Wang et al. 2007; Younger et al. 
2007). 

Soon after the first observations with SCUBA, surveys were 
being conducted at 1200 with MAMBO (Bertoldi et al. 2000) 
resulting in samples of sources that were thought to be identical to 
SMGs, except selected at mm wavelengths. Compared to submm 
surveys, however, observations at mm wavelengths should be sen- 
sitive to far-IR-luminous galaxies out to higher redshifts, due to 
the fc-correction, as well as to systems with cooler dust tempera- 
tures. As a result, at z > 3 the 850-/1200-/im flux ratio becomes 
a relatively strong function of redshift for a typical starburst far- 
infrared/mm spectral energy distribution (SED), and may therefore 
(if the SED is known) be used as a crude redshift estimator or - if 
the redshift is known - as an indicator of the shape (i.e. temperature 
and spectral index) of the far-infrared/mm SED. 

This fact was exploited by Eales et al. (2003; hereafter E03), 
who obtained 850-/im photometric measurements of a sample of 
1200-/im-selected sources and found a significant fraction to have 
very low 850/1200-/im flux ratios (S's50Mm/S'i200nm £ 2). They 
argued that this was either indicative of SMGs at z ^> 3 or - 
since a large fraction of their sample was radio-identified with 
radio-to-submm flux ratios consistent with z ^ 3 - SMGs having 
fundamentally different dust emission properties than local galax- 
ies. The existence of a significant number of these SDOs, which 
we define here as sources with <S , 85o M m/Si200/jm ^ 2, was ques- 
tioned (albeit not ruled out) by Greve et al. (2004; hereafter G04), 
who performed the first unbiased comparison of 850- and 1200- 
^im-selected sources, using catalogues extracted from SCUBA and 
MAMBO maps. The degree of overlap between source popula- 
tions extracted from submm and mm surveys (and in particu- 
lar the idea that mm observations are tracing a significant sub- 
population of more distant and/or cooler sources) remains contro- 
versial, therefore, and is further complicated by the difficulty of 
comparing (sub)mm surveys with different depths and noise prop- 
erties. Nonetheless, settling this issue is important; regardless of 
the outcome, it tells us about the population of dusty galaxies at the 
highest redshifts, with significant implications for our understand- 
ing of galaxy formation and evolution. 

In this paper we present an unbiased 1200-/im MAMBO sur- 
vey of the northern field of the Great Observatories Origins Deep 
Survey (GOODS-N - Giavalisco et al. 2004). As well as presenting 
the 1200-/im map and source catalogue, we also perform a rigor- 
ous comparison between our data and the existing 850-/iin SCUBA 
map of this region (Borys et al. 2003; Pope et al. 2005), in order to 
address the issues outlined above. In a follow-up paper, we shall 
exploit other high-quality data in the GOODS-N field to explore 
the multi-wavelength properties of 1200-/im selected sources. 

Throughout this paper we have adopted a flat cosmology with 
Q,m = 0.27, fU = 0.73, and H = 71 kms" 1 Mpc -1 (Spergel et 
al. 2003). 



2 OBSERVATIONS AND DATA REDUCTION 

A 1200-/im (250-GHz) map covering the entire GOODS-N region 
was obtained using the 117-channel bolometer array, MAMBO 
(Rreysa et al. 1998), on the IRAM 30-m telescope on Pico Veleta in 
Spain. The data were obtained during the winter semester 2005-06 
in good weather conditions, meaning 250-GHz opacities of <0.3 
and stable atmospheric conditions (low sky-noise). 

In order to map the entire field as uniformly as possible, a 
number of regular grid positions (2 arcmin apart) were observed. 
Specifically, each such grid position was observed twice, using 
a standard on-the-fly 300arcsec x 320arcsec Az-Alt scan-map. 
With a scanning speed of 5 arcsec s _1 in Az and an elevation spac- 
ing between each sub-scan of 8 arcsec, a scan-map takes ~43 min 
to complete, allowing for a 3-s turnaround after each sub-scan. To 
avoid any systematics and to minimise the effects of sky rotation, 
we tried to observe each grid position once when the field was ris- 
ing and once when it was setting. Different wobbler throws (36- 
45 arcsec) and scan directions (increasing or decreasing elevation) 
were also employed for different scan-maps. In total, 53 scan-maps 
were obtained. 

After each scan-map, a pointing observation on a nearby, 
bright point source was carried out and, in order to be extra cau- 
tious, a focus correction was usually applied after each scan (some- 
times after two scans if the focus appeared stable). Since the 
weather conditions were generally stable, it was sufficient to per- 
form a skydip every 2hr in order to monitor the 250-GHz atmo- 
spheric opacity. Skydips were carried out at the same azimuth as the 
preceeding grid-point observation, ensuring that the opacity was 
measured for the same part of the sky as the science observations. 

The data were reduced with the MOPSIC software - an im- 
proved version of the original MOPSI package (Zylka 1998) - using 
standard pipeline scripts designed to deal with MAMBO data. This 
involved many of the standard steps in reducing bolometer data, 
i.e. de-spiking, flat-fielding, removal of correlated noise across the 
array (see G04 for details), and correction for atmospheric opacity 
by interpolating between the skydip measurements. The timestream 
data were rebinned onto a 1 arcsec x 1 arcsec grid using a shift-and- 
add technique, which results in each positive source being brack- 
eted by two negative counter-images, each with half of the peak 
positive intensity, located one wobbler throw away, i.e. the standard 
'double-difference' or 'triple-beam' pattern. However, the differ- 
ent wobbler throws and scanning directions employed for different 
scans ensure that the negative images are averaged out, resulting in 
a very clean image (although with some residual effects, discussed 
in §[ID. 

The final 1200-/im signal-to-noise (hereafter S/N) map and 
corresponding noise map are shown in Fig. Q] The area covered 
is 287 arcmin 2 . The map exhibits a fairly uniform r.m.s. noise of 
~0.7 mJy beam - over most of its area (apart from at the edges 
where the noise is larger, ^0.8 mJy beam" 1 ). 

3 THE GOODS-N 1200-MICRON SOURCE CATALOGUE 
3.1 Source extraction 

In order to extract sources, a point spread function (PSF) was fitted 
to each pixel in the map and a S/N ratio was ascribed, based on 
the local noise. This is similar to the noise- weighted convolution 
technique applied to MAMBO maps by G04 and to SCUBA maps 
by Borys et al. (2003) and Coppin et al. (2006). Unlike SCUBA, 
where one can configure the nod/chop such that it remains constant 
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Figure 1. Top: The 1200-/xm MAMBO S/N map of the GOODS-N region. Sources extracted at JS3.5 <r (prior to deboosting, see i|3.2.3t are indicated as white 
circles. Bottom: The corresponding r.m.s. noise map. The thin black contours correspond to a = 0.7, ... , 1 . 2 mjy beam - 1 . The thick black outline shows the 
full extent of the GOODS-N field as observed with Spitzer/YRAC and the thick dashed line shows the extent of the SCUBA map (Borys et al. 2003; Pope et al. 
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GOODS-N CENTRE: R.A. '2:37:00.0 DEC. +62:'3:00.0 
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Figure 2. a) Variations in the PSF at a collection of grid points across the map. The negative parts of the PSF (shown in white) are caused by the wobbler 
throw used in the on-the-fly mapping mode on the IRAM 30-m telescope. Due to sky rotation and variations in the chop throw (between 36 and 45 arcsec) 
the negative residuals are smeared out in a ring around the main beam, b) Azimuthally averaged PSF profiles. The gray-shaded region shows the radially 
averaged profiles of all the PSFs plotted on the left, and the dotted line is the average of those profiles. The solid line is a best-fit Gaussian with a fixed FWHM 
of 11.1 arcsec. 



on the sky, MAMBO experiences sky rotation; as a result, the PSF 
varies across the map. 

We have gone to great lengths to average out the negative chop 
patterns (see §[2](, but we still expect to see some effects of the chop 
throw in the map. To quantify this, we created fake, noiseless time 
streams of data, containing a number of point sources at certain po- 
sitions in the map, and then subjected these data to exactly the same 
data-reduction pipeline as the real data, using the astrometrical in- 
formation (coordinates, chop throws etc.) of the real data to create 
a map which shows how the PSF varies across the map (Fig.|2^). 

Ideally, one would want to use the PSF (including its negative 
parts) as the filtering kernel in the source-extraction procedure, as 
has been done with SCUBA maps (e.g. Scott et al. 2002, hereafter 
S02; Borys et al. 2003; Coppin et al. 2006). In our case, that would 
entail computing a separate PSF for every pixel in the map, or at 
least some sub-region of the map where the PSF is approximately 
constant. However, as seen in Fig. [2^, the negative features of the 
net PSF are fairly low contrast, spread out over a large area, making 
it difficult (and also largely unnecessary) to use them explicitly in 
the source-extraction procedure. 

This also suggests that the optimal method of combining the 
data, using an explicit PSF for each individual time sample (see 
Borys et al. 2003), would produce only marginally better results, 
at the expense of enormously increased complexity. Indeed, we 
used a few sources to show that the improvement was modest. We 
therefore decided to use a simple Gaussian (FWHM, 11.1 arcsec) 
as an approximation to the PSF. Although the 30-m telescope's 
beam is formally 10.7 arcsec FWHM at 1200 /im, we have adopted 
a slighter broader PSF in order to account for the typical pointing 
error (~3 arcsec r.m.s.). That this is, in fact, a good approxima- 
tion is demonstrated in Fig. |2p where we have plotted the mean 



azimuthally averaged PSF profiles, along with a Gaussian profile 
with FWHM =11.1 arcsec. While half of the signal is contained in 
the negative residuals, it is seen that the depth of the negative holes 
reach only about 5 per cent of the primary PSF height. This proce- 
dure yielded 20 sources with S/N ^ 4 and a further 27 sources with 
3.5 < S/N < 4 (see TableQ}. 

3.2 Monte Carlo simulations 

3.2.1 Spurious sources 

In order to assess the degree to which the extracted source cata- 
logue is contaminated by spurious noise peaks, we randomised the 
bolometer positions on the array in order to create a scrambled map 
with no sources, with noise properties similar to that of the real 
map. In total, 500 such scrambled maps were produced, each pro- 
cessed by the data-reduction pipeline in exactly the same way as the 
real data. Fig. [3] shows the average number of spurious detections 
in the scrambled maps, when subjected to our source-extraction 
technique. The simulations suggest that for ^4-ct sources extracted 
from the 287 arcmin 2 surveyed, less than two are expected to be 
spurious, while we can expect about seven spurious sources in 
the 3.5-ct catalogue. This is in agreement with similar simulations 
by G04. Fig. [3] also shows how the expected number of spurious 
sources is larger than for purely uniform Gaussian noise, confirm- 
ing that it is important to carry out the simulations with realistic 
noise. 

A further check on the source robustness can be made by ex- 
tracting sources from an inverted map, in other words, carrying out 
the same extraction procedure for negative sources, which of course 
can only be real when created by freak alignments of sources and 
chop throws. In total, we find 33 sources above 3.5 a in the inverted 
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Figure 3. Diamonds show the average number of (positive) spurious de- 
tections as a function of S/N ratio, based on simulated maps with realis- 
tic noise properties. The dashed line shows the expected number of spuri- 
ous sources, assuming purely Gaussian noise and a Gaussian beam with a 
FWHM of 11.1 arcsec. 



map. All but six are associated with the negative beams of bright 
positive sources, consistent with expectations. 



3.2.2 Completeness and positional offsets 

Catalogue completeness varies as a function of the source flux, be- 
coming smaller at lower flux values. In addition, random noise fluc- 
tuations result in the extracted positions being offset from the true 
source positions; again, the effect becomes increasingly severe at 
lower fluxes. The effect of the noise on the completeness of the sur- 
vey and on the positional offsets was assessed using Monte Carlo 
simulations. Artificial sources with fluxes in the range 1-12 mjy 
(in steps of 0.25 mjy) were added to the science map, one at a time, 
at random positions, and subsequently extracted using exactly the 
same technique used in § 13.11 This procedure was repeated 50 times 
for every flux bin, each time at a random position in the map. Each 
injected source was a scaled version of the local PSF and as a result 
the shape and depth of the negative lobes around the sources was 
a function of the map position. We also excluded artificial sources 
which were within half a beam width of any real source, hence this 
definition of completeness is intended only to assess the effects of 
noise, without including the effects of confusion (Blain, Ivison & 
Smail 1998). 

The resulting completeness curve and positional offsets are 
shown in Fig. [4] The simulations suggest that the survey is close 
to complete down to a flux density of about 3 mjy, and drops pre- 
cipitously below 2.5 mjy. The lower panel of Fig.|4]shows that the 
positions of sources with flux densities larger than 2.5 mjy are de- 
termined to within 3 arcsec. Of course, in addition to the positional 
offsets shown in Fig. [4] the data will contain offsets due to pointing 
errors. The completeness fractions and positional offsets we obtain 
here are very similar to those found by G04. 




2 4 6 8 10 12 
S in [mjy] 

Figure 4. Top: The completeness function of our survey, derived as the per- 
centage of input sources which are recovered by the source extraction as a 
function of the source flux density (but not including the effects of confu- 
sion). The solid curve represents a fit to the phenomenological expression: 
/(Sj n ) = 1 - exp(A(5 in - B) c ). Bottom: The positional offsets be- 
tween the actual positions of the input sources and their extracted position. 
A fit to the function /(Sj n ) = 1 + Aexp(B SP) is shown by the solid 
curve. 



3.2.3 Flux deboosting 

For any flux-limited survey in the low-S/N regime (as is the case 
for all submm/mm surveys undertaken to date), we must worry 
about the superposition of sources with positive noise peaks in the 
map: 'flux boosting' (S02). This effect is particular significant in 
the submm/mm waveband where the source number counts have a 
steep negative slope, so that the noise excursions above a threshold 
do not balance excursions below the threshold. 

We chose to follow the Bayesian flux deboosting approach de- 
veloped by Coppin et al. (2005) to deal with SCUBA data. The phi- 
losophy of this method is to assume a priori knowledge of the pixel 
brightness distribution, N(S P ), in the noiseless map. If the noise in 
the map is Gaussian, the probability of measuring a flux density 
S m ± c m at some pixel in the map when its true flux is S p , is 
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Figure 5. Left: The prior flux density probability distribution, N(S P ), of 1,000 simulated maps. Whilst strongly peaked around zero, the distribution is skewed 
towards positive fluxes due to the presence of sources in the maps. The flux density distribution of the real map is shown in grey. Middle: Observed flux density 
probability distribution (dotted curve) assuming Gaussian errors, Pm(S p ), and the posterior flux density probability distribution (dark solid curve), P(S P ), 
for the source MMJ12371 1+621328 (GN 1200.3), which is well detected. Right: An example of a source, MM J123541+621150 (GN 1200.35), which is 
less well detected. In this case the posterior flux density distribution peaks at zero and the posterior probability for having zero flux density is larger than 5 per 
cent. Hence this source does not meet our detection criterion, after deboosting. 



Pm(S p ) OC exp 



2ai 



(1) 



We can then use Bayes' theorem to derive the probability that the 
true flux of the pixel is S p : 



P(S P ) = 



N(S P )P M (S P ) 
P(M) 



(2) 



where P(M) is simply the normalisation constant. The a priori 
flux density distribution, N(S P ), was generated by making noise- 
less maps onto which sources selected from a model of the num- 
ber counts were randomly distributed, and then sampled with the 
adopted chop pattern. In our case, the chop pattern varies across 
the map so the overall N(S P ) distribution was derived from noise- 
less maps created with exactly the same scanning strategy as the 
real map. For the number counts, we adopted a Schechter function 
oftheformdAVdS' oc (S/S )~ a exp(-S/S ), where So = 4mJy 
and a = 2.3; this provides a good fit to the 1200-^im number 
counts by G04. The N(S P ) distribution derived from 1,000 such 
maps is shown in the left panel of Fig. [5] We checked that the use 
of moderately different number count models has little effect on 
our results; one has to considerably alter the faint-end slope in or- 
der to make significant changes. Examples of deboosting are given 
in Fig. [5] which shows the probability distribution of the measured 
flux Pm(S p ) (assuming it is Gaussian) and the corresponding de- 
boosted flux distribution P(S P ) for the two sources, GN 1200.3 and 
GN 1200.35. 

Following Coppin et al. (2006), only sources with a posterior 
probability less than 5 per cent of having a flux density less than 
or equal to zero were retained in the final source catalogue. This 
means that of the two examples just employed, GN 1200.3 makes 
it into the final catalogue, whereas GN 1200.35 does not. Adopt- 
ing this criterion leads to a total of 17 sources being discarded, all 
of them originally detected at below 4 a. The final revised source 
list is indicated through column seven of Table Q] We adopt the 
deboosted flux density as the value at which the P(S P ) distribu- 
tion peaks, with the errors calculated as the 68 per cent confidence 
bounds on the non-Gaussian deboosted flux distribution. On aver- 
age, the deboosting reduces the observed fluxes by about 0.7 mjy, 
as shown in Fig. [6] In addition, the deboosting results in larger 
photometric errors and non-Gaussian probability distributions (see 
Fig. O, particularly in cases where the photometric errors are al- 
ready large. 



The final catalogue of deboosted MAMBO sources 

consists of 30 sources. In ^3.2.11 we found that seven of the 
47 ^3.5-<t sources recovered prior to deboosting were spurious. 
Hence, of the sources in the deboosted list we expect ~7 x (30/47), 
or around four, to be spurious. This is almost certainly an upper 
limit since deboosting is designed to preferentially remove spuri- 
ous sources. In Table[T]we list all of the original ^3.5-cr sources, 
because many of those which did not pass our deboosting criteria 
still have a reasonably high chance of being genuine and hence may 
be useful for follow-up studies at other wavelengths. 



4 COMPARISON WITH THE SCUBA GOODS-N SURVEY 

The SCUBA 850-^m survey of the Hubble Deep Field North 
(HDF-N), which lies at the centre of the GOODS-N field, was 
one of the first extragalactic fields to be observed at submm wave- 
lengths (Hughes et al. 1998; Serjeant et al. 2003) and to this day re- 
mains one of the deepest blank-field submm maps. A much larger, 
shallower 850-/^m image, covering most of the GOODS-N field, 
was presented by Borys et al. (2002) and this has been supple- 
mented by several SCUBA jiggle maps (Holland et al. 1999) of 
other parts of GOODS-N (Borys et al. 2003, 2004; Wang et al. 
2004; Pope et al. 2005; Pope et al. 2006). These data were com- 
bined into the so-called 'super-map' of GOODS-N by Borys et al. 
(2003), then updated in a series of subsequent papers. The resulting 
map combined the available data optimally, with the purpose of ex- 
tracting point sources. The final version of the 850-/im source cat- 
alogue was presented by Pope et al. (2006), with three new sources 
(GN 850.39, GN 850.40 and GN 850.41) presented by Wall, Pope 
& Scott (2008). 

The AzTEC camera (Wilson et al. 2008) was used to survey 
the GOODS-N region at 1 100-/im (Perera et al. 2008). A compari- 
son of the MAMBO and AzTEC millimetre maps will be explored 
in a future paper (Pope et al. in preparation). 

In this section we will compare the 850-/im super-map with 
our new 1200-/im image, with two goals: first, to extract a robust 
sample of SMGs detected jointly at 850 and 1200 /im; second, to 
compare the degree of overlap between the populations selected at 
each wavelength. 

The optimal method for extracting sources of arbitrary spec- 
tral shape is to extract sources from the combined S/N map, thereby 
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Table 1. 1200-/im MAMBO sources extracted from the GOODS-N field at a significance ^4<r, with a supplementary list of sources with 3.5 ^ S/N < 4 
(italics). Sources which pass our deboosting criterion are assigned a deboosted flux density in the 7th column. These 30 sources make up our robust, deboosted 
catalogue, of which we expect no more than 3—4 to be spurious ( i]3.2.3t . The three last columns list 1.4-GHz flux densities of all robust radio sources within a 6- 
arcsec radius of each MAMBO source, their mm-radio positional offset and the statistical significance of each mm-radio association, respectively. Associations 
with P < 0.05 are considered robust and are listed in boldface (see also § 15.41 . 
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GN name 


R.A. (J2000) 


Dec. (J2000) 




11 -t cr 1200^m 

/mjy 


S/N 


Sl200Aim (deboosted) 
/mjy 


5l.4GHz 
/ply 


offset 
/arcsec 
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MM J lz_1 / 1 1 -\-oZZZ 1 1 


uJN tzUU.t 
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9 


3 + 0.9 




9.3 i 0.5 


58 + 16 






MMJ123633+621407 


GN 1200.2 


12:36:33.2 


+ 62:14:07 


4 


9 ± 0.7 


6.95 


4.4 ± 0.7 


34 ± 4 


2.28 


0.017 


MM J12371 1+621328 


GN 1200.3 


12:37:11.2 


+62:13:28 


3 


9 ± 0.6 


6.50 


3.6 ± 0.6 


149 ± 14 
62 + 6 


3.16 
5.92 


0.009 
0.047 


MMJ123730+621258 


GN 1200.4 


12:37:30.8 


+62:12:59 


4 


2 ± 0.7 


6.04 


3.7 ± 0.7 


114 ± 7 


0.15 


0.00006 


MMJ123626+620604 


GN 1200.5 


12:36:26.9 


+62:06:04 


4 


7 ± 0.9 


5.16 


3.8 ± 1.0 


46 + 8 


3. IS 


0.025 


MMJ123714+621827 


GN 1200.6 


12:37:14.1 


+62:18:28 


3 


6 ± 0.7 


5.08 


2.9 ± 0.8 


627 ± 8 


2.33 


0.002 


MM J 1 23607+620855 


GN 1200.7 


12:36:07.8 


+62:08:56 


3 


4 ± 0.7 


4.XI) 


2.8 ± 0.8 








MMJ123702+621948 


GN 1200.8 


12:37:02.9 


+62:19:48 


3 


3 ± 0.7 


4.67 


2.7 ± 0.8 








MMJ123644+621940 


GN 1200.9 


12:36:44.6 


+ 62:19:40 


3 


3 ± 0.7 


4.65 


2.7 ± 0.8 








MMJ123737+621456 


GN 1200.10 


12:37:37.8 


+62:14:56 


2 


8 ± 0.6 


4.61) 


2.2 ± 0.7 








MMJ123609+620648 


GN 1200.11 


12:36:09.7 


+62:06:48 


3 


6 ± 0.8 


4.48 


2.8 ± 0.9 








MMJ123550+621041 


GN 1200.12 


12:35:50.3 


+ 62:10:42 


4 


5 ± 1.0 


4.48 


3.5 ± 1.2 








MMJ123631+621715 


GN 1200.13 


12:36:31.8 


+ 62:17:15 


3 


± 0.7 


4.25 


2.2 ± 0.8 








MMJ123635+620714 


GN 1200.14 


12:36:35.6 


+ 62:07:14 


3 


± 0.7 


4.25 


2.2 ± 0.8 








MMJ123554+621338 


GN 1200.15 


12:35:54.0 


+62:13:38 


3 


± 0.7 


4.21 


2.2 ± 0.9 


55 + 8 


5.13 


0.043 


MMJ123723+622049 


GN 1200.16 


12:37:23.7 


+62:20:49 


3 


3 ± 0.8 


4.15 


2.3 ± 1.0 








MMJ123736+621103 


GN 1200.17 


12:37:36.0 


+ 62:11:04 


2 


9 ± 0.7 


4.12 


2.1 ± 0.9 


42 + 7 


5.89 


0.057 


MMJ123556+621726 


GN 1200.18 


12:35:56.0 


+62:17:26 


3 


3 ± 0.8 


4.07 


2.3 ± 1.0 








MMJ123645+620823 


GN 1200.19 


12:36:45.1 


+62:08:24 


2 


8 ± 0.7 


4.02 


2.1 ± 0.9 








MMJ123749+621558 


GN 1200.20 


12:37:49.3 


+62:15:59 


2 


8 ± 0.7 


4.01) 


2.1 ± 0.9 








3 . 5 ^ 0" < 4 detections 


MM J 123720+621 106 


GN1200.21 


12:37:20.8 


+62:11:06 


2 


8 + 0.7 


3.95 


'2.0 + 0.9 








MM J12371 1+620652 


GN 1200.22 


12:37:11.2 


+62:06:52 


3 


9+1. 


3.92 










MMJ123717+620800 


GN1200.23 


12:37:17.4 


+62:08:00 


3 


1 +0.8 


3.91 


2.2 + 1.1 


97 ± 19 


4.24 


0.021 


MMJ123755+621057 


GN 1200.24 


12:37:55.6 


+62:10:58 


3 


5 + 0. 9 


3.91 










MMJ123650+621217 


GN1200.25 


12:36:50.9 


+62:12:17 


2 


7+0.7 


3.88 


2.0 + 0.9 








MM J 123643 +622350 


GN 1200.26 


12:36:43.9 


+62:23:51 


4 


7+1.2 


3.87 










MM 1123607+621137 


GN1200.27 


12:36:07.2 


+62:11:37 


2 


7+0.7 


3.82 


2.0+1.0 








MMJ123605+621650 


GN1200.28 


12:36:05.8 


+62:16:50 


2 


7+0.7 


3.80 


2.0 + 1.0 








MMJ123711+62121S 


GN1200.29 


12:37:11.8 


+62:12:18 


2 


7+0.7 


3.80 


2.0+1.0 








MM J123638+622049 


GN1200.30 


12:36:38.8 


+62:20:50 


3 


+ 0.8 


3.78 










MM 1123803+621446 


GN1200.31 


12:38:03.4 


+62:14:46 


3 


+ 0.8 


3.77 










MM J 1237 11 +62 1924 


GN1200.32 


12:37:11.8 


+62:19:24 


2 


6 + 0.7 


3.75 


1.9 + 1.0 








MMJ123641+621129 


GN1200.33 


12:36:41.9 


+62:11:30 


2 


6 ± 0.7 


3.74 


1.9 + 1.0 








MM 1123752+621428 


GN1200.34 


12:37:52.0 


+62:14:28 


2 


6 ± 0.7 


3.70 










MM 1123541+621150 


GN1200.35 


12:35:41.7 


+62:11:50 


4 


1+1.1 


3.69 










MM 1123802+621940 


GN1200.36 


12:38:02.8 


+62:19:40 


4 


4 + 1.2 


3.68 










MM 1123656+621647 


GN1200.37 


12:36:56.8 


+62:16:47 


2 


6 + 0.7 


3.65 










MM 1123540+621246 


GN1200.38 


12:35:40.3 


+62:12:47 


3 


7+1.0 


3.65 










MM 1123714+621842 


GN1200.39 


12:37:14.1 


+62:18:43 


2 


6 + 0.7 


3.64 










MM 1123745+621 129 


GN1200.40 


12:37:45.4 


+62:11:30 


2 


5 ± 0.7 


3.61 










MM 1123621+621705 


GN1200.41 


12:36:21.5 


+62:17:06 


2 


5 ± 0.7 


3.55 




138 ± 7 
50 ± 7 


2.S6 
5.24 


0.008 
0.046 


MM 1123700+621453 


GN1200.42 


12:37:00.6 


+62:14:54 


2 


1 + 0.6 


3.55 


1.4 + 0.8 


53 ± 13 


5.84 


0.051 


MM 1123712+621258 


GN1200.43 


12:37:12.7 


+62:12:59 


2 


1 ± 0.6 


3.55 


1.4 + 0.8 








MM 1123659+621619 


GN1200.44 


12:36:59.2 


+62:16:19 


2 


5 ± 0.7 


3.54 










MM 1123653+622401 


GN1200.45 


12:36:53.4 


+62:24:01 


4 


2 + 1.2 


3.53 










MM 1123655+621408 


GN1200.46 


12:36:55.5 


+62:14:09 


2 


1 + 0.6 


3.51 










MM 1123557+621250 


GN 1200.47 


12:35:57.6 


+62:12:50 


2 


5 + 0.7 


3.51 











using all the information from the 1200- and 850-^im data. How- 
ever, some care must be exercised, for at least two reasons. First, 
one is digging into the confusion noise in both images and hence 
non-Gaussianity or systematic artefacts in the maps may start to 
have effects that were not previously evident. Second, if one had 
strong prior information about the expected Sssonm/S'wooMm ratio 
for the emitters then one would adopt a different weighting scheme. 

The source-extraction method applied to the two maps indi- 
vidually (see Borys et al. 2003) is optimal, with the caveats given 
in § 13.11 for finding point sources, under the assumption that the 
noise distributions in both maps are Gaussian and uncorrelated. 



The SCUBA super-map can be thought of as the answer to the 
question: 'what is the best estimate of the flux density of a point 
source centred on each pixel?'. The super-map combines data taken 
in different observing modes, with different chop-throws, etc., and 
also takes into account the beam convolution in the process of fit- 
ting to the PSF. It follows that an efficient way of extracting point 
sources, using the information in both maps, is simply to treat the 
MAMBO image as additional data. These can then be combined 
with the SCUBA data to provide an overall S/N 'map', which gives 
the combined S/N level for a point source centred on each pixel. 
The relevant expression for each pixel is: 
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CO 



2.0 



1.5 



1.0 



CO 



0.5 



0.0 



• SNR>4.0a 

° 4.0a >SNR> 3.5a 



2 3 4 5 6 

S ra [mjy] 

Figure 6. The difference between the raw and the deboosted 1200-/xm 
flux densities as a function of the raw flux densities for sources in the fi- 
nal deboosted catalogue. On average, the raw flux densities are reduced by 
~0.7 mJy as a result of the deboosting. 



S = 



(3) 



where ctm and as are the MAMBO and SCUBA noise maps, re- 
spectively. The corresponding noise map is 



M + °S- 



(4) 



For this method to work straightforwardly, the two maps need to be 
astrometrically aligned and to have the same pixel scale and reso- 
lution. To this end, both maps were binned onto a 1-arcsec 2 pixel 
grid and we checked that no shift between the grids was required 
in order to maximise the correlation, then the MAMBO image was 
smoothed to the SCUBA resolution. We stress that this method only 
improves on the robustness of our detections, it does not provide es- 
timates of the 850- and 1200-/im flux densities, which have to be 
measured from the separate maps. 

Fig.|7]shows the combined 850- and 1200-/iin S/N and r.m.s. 
noise maps obtained in this way, and Table [2] lists the sources ex- 
tracted from the combined map. 

A total of 13 (33) sources are detected above 5 a (4 a) in 
the combined map. In the bottom panel of Fig. [7] and in Table [2] 
we compare our dual-map source list with the SCUBA-only and 
MAMBO-only source lists. Any source from the combined map 
which has a >3.5-a deboosted MAMBO (SCUBA) source within 
lOarcsec (V8 2 + 6 2 arcsec, where 8 and 6 arcsec are the 95-per- 
cent confidence limit for SCUBA and MAMBO sources (Ivison et 
al. 2005; Pope et al. 2006) is deemed to have been associated at 
1200 ^trn (850 ^m). 

How many sources have robust detections at both 850 and 
1200 ^im? Only six sources from the combined list, four of which 
are detected at ^5 a in the combined map, and one at 4 a. The sixth 
source, GN 1200.42 (GN 850.32), is detected at ~3.8a. In cases 
with no detection at 850 or 1200 /jm, we have determined 3-a up- 
per flux density limits (Table |2). Finally, we note that one source 
detected at >4 a in the combined map, MM J12371 1+621245, was 
not previously known. 

How big a fraction of the MAMBO and SCUBA sources are 
recovered in the combined map? All five (14) of the ^5-a (^4-a) 
MAMBO sources which lie within the SCUBA region are recov- 
ered in the combined map. Conversely, all of the ^5-a SCUBA 



sources are detected in the combined map, while for the ^4-a 
SCUBA sources the corresponding fraction is 18/24 (75 per cent). 

If we instead use Table [2] to answer the question 'how big a 
fraction of the MAMBO sources have robust associations in the 
SCUBA map?', we find that apart from the four ^5-a MAMBO 
sources with SCUBA counterparts (a 100-per-cent recovery rate 
since the fifth ^5-a MAMBO source, GN 1200.5, lies outside the 
SCUBA region) none of the other ^4-a MAMBO sources have 
been recovered at 850 /im. This is clearly at odds with the high frac- 
tion of associations found for ^4-a MAMBO sources by G04 in 
ELAIS N2 and Lockman Hole (~30-50 per cent) and the similarly 
high rate of associations for bright 1.1 -mm BOLOCAM sources 
in the Lockman Hole (~75 per cent - Laurent et al. 2006). While 
positional uncertainties may go some way towards explaining the 
low fraction of associations in GOODS-N, the dominant factor is 
probably the mis-match in sensitivity of the MAMBO and SCUBA 
surveys in this region. 

Finally, we note that GN 1200.4 - one of the brightest 
(<Si200Mm — 4.2 ± 0.7 mjy) and most significant (>5 a) sources 
in our sample - is not detected at 850 /im and is thus a submm 
drop-out (SDO) candidate, although it does lie in a relatively noisy 
part of the SCUBA map. 



5 THE 850-/1200-MICRON FLUX DENSITY RATIO 

In this section we shall examine the 850-/im to 1200-/^m flux 
density (SsBo^m / Si200^m) distribution of the (sub)mm sources in 
GOODS-N. 



5.1 Sources detected by both SCUBA and MAMBO 

We begin by considering the sources which are robustly detected by 
both MAMBO and SCUBA, i.e. sources which are in both of the 
deboosted source catalogues and detected at ^4 a in the combined 
map. In total, only five such sources were found (see Table |2) and 
their resulting SssoMm/S^oo/jm versus 850-/im flux density distri- 
bution is shown in Fig. [8] 

For comparison Fig. [8] also shows the robust 
MAMBO/SCUBA sample from E03, which are based on 
SCUBA photometry of MAMBO sources extracted from surveys 
of the Lockman Hole, the NTT Deep Field and the Abell2125 
cluster field (Bertoldi et al. 2000; Dannerbauer et al. 2002; 2004). 
Values of SssoAim and 5i200nm m E03 did not require deboosting 
as they were measured using 'photometry mode' (see Ivison et 
al. 1998b) at radio and/or mm interferometric positions. We also 
compare with the MAMBO/SCUBA samples derived from blank- 
field surveys of the ELAIS N2 and Lockman Hole fields by S02 
and G04. For the five MAMBO/SCUBA associations found in the 
ELAIS N2 field, we derived Si200nm using the 'deboosting curve' 
given in G04, while a uniform deboosting factor of 15 per cent was 
adopted for 5g5o Mm (S02). In both cases, the deboosting factors 
were within the total flux calibration errors. Using a dual-survey 
extraction technique, Ivison et al. (2005 - hereafter 105) identified 
eight MAMBO/SCUBA associations in the Lockman Hole - in 
addition to the eight already presented by G04 - and we have 
included those in our comparison here. For all 16 associations in 
the Lockman Hole we have adopted Ssso^m (deboosted) from the 
SHADES survey (Coppin et al. 2006) which has a larger amount 
of data on the Lockman Hole than the UK 8-mJy Survey (S02). 
Deboosted Si2oo^m came from G04. 

Although we have only identified five ^4-a 
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GOODS-N CENTRE: R.A. 12:37:00.0 DEC. +62:13:00.0 




500 -500 

Acx (arcsec) 

GOODS-N CENTRE: R.A. 12:37:00.0 DEC. +62:13:00.0 




Acx (arcsec) 

Figure 7. S/N image (top) and r.m.s. noise map (bottom) resulting from combining the 850-/xm SCUBA map (Borys et al. 2003; Pope et al. 2005) and 1200-/^m 
MAMBO map (this work) according to eqs|5]and|4] White circles indicate sources with significance from the combined image. For the sake of clarity, 

these sources have not been numbered in the bottom panel. Green triangles indicate the deboosted ^3.5-cr MAMBO sources in Table[T] yellow crosses are the 
deboosted ^3.5-cr SCUBA sources, as given by Pope et al. (2006), with the addition of GN 850.39, .40 and .41 from Wall, Pope & Scott (2008). 
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Table 2. Sources extracted at ^4 cr from the combined 1200- and 850-/im image (eqs[3]and[4), with a supplementary list of sources extracted at 3.5 ^ S/N < 4 
(italics). The indices in parentheses refer to the numbering of the MAMBO catalogue in Table 1, except when there is no MAMBO association in which case 
we have started the index at 100. The same numbering is used in Fig. [7] The MAMBO and SCUBA associations in the second column from the left are those 
given in Table[T]and in Pope et al. (2005), respectively. Deboosted 1200- and 850-/xm flux densities are given in parentheses in columns 6 and 7. In cases where 
no robust associations have been identified, we list the 3-cr upper flux density limit. Columns 8, 9 and 10 list Si.4GHz f° r a U robust radio identifications within 
8 arcsec of each MAMBO/SCUBA association in the combined map, the radial offset and the statistical significance of each radio identification, respectively. 
Identifications with P < 0.05 are considered robust and are given in boldface. 



Name 


GN name 


R.A. (J2000) 


Dec. (J2000) 


S/N 


Sl200/im 1 <T 1200Mm 




Si 4GHz 


offset 


P 
















/mjy 




/mjy 


/li}y 






(100)MMJ123652+621226 


GN 850.14 


12:36:52.2 


+62 


12:26 


14.7 




^ 2.1 


5 


9 ± 0.3(5.9) 


< 16 






(1) MM J12371 1+622212 


GN1200.1/GN 85020 


12:37:11.5 


+62 


22:12 


13.5 


9 


3 ± 0.9 (9.3) 


20 


3 ± 2.1 (20.3) 


58 ± 16 


2.63 


0.016 


(101)MMJ123656+621205 


GN 850.15 


12:36:56.4 


+62 


12:05 


8.75 




^ 2.0 


3 


7 + 0.4(3.7) 


42 + 6 


3.02 


0.024 


(2) MM 1123633+621407 


GN1200.2/GN 850.10 


12:36:33.2 


+62 


14:07 


7.85 


4 


9 ± 0.7(4.4) 


11 


3 ± 1.6(11.3) 


34 + 4 


2.28 


0.017 


(3)MMI123711+621328 


GN 1200.3/GN 850.39 


12:37:11.1 


+62 


13:28 


7.16 


3 


9 ± 0.6 (3.6) 


7 


4 ± 2.0 (5.2) 


149 ± 14 
62 + 6 


3.42 
6.57 


0.010 

0.053 


(4)MMJ123730+621300 


GN 1200.4 


12:37:30.7 


+62 


13:00 


6.14 


4 


2 ± 0.7(3.7) 




^ 15.6 


114 ± 7 


1.23 


0.002 


(6)MMI123714+621827 


GN1200.6/GN 850.40 


12:37:14.0 


+62 


18:27 


(,.()() 


3 


6 ± 0.7(2.9) 


13 


1 ± 2.7(10.7) 


627 ± 7 


1.16 


0.0005 


(102)MMJ123701+621148 


GN 850.17 


12:37:01.1 


+62 


11:48 


6.00 




^ 2.0 


3 


9 ± 0.7(3.9) 


102 ± 12 


3.31 


0.014 


(103)MMJ123621+621707 


GN 850.7 


12:36:21.5 


+62 


17:07 


5.83 




$ 2.0 


8 


9 ± 1.5(8.9) 


138 ± 7 


2.09 
4.59 


0.005 
0.039 


(104)MMJ123637+621158 


GN850.11 


12:36:37.1 


+62 


11:58 


5.83 




^ 2.0 


7 


± 0.9(7.0) 


26 + 5 


3.19 


0.028 


(105)MMJ123650+621315 


GN850.13 


12:36:50.7 


+62 


13:15 


5.25 




^ 2.1 


1 


9 ± 0.4(1.9) 


47 + 6 


7.27 


0.068 


(106) MM J1236 18+62 1551 


GN 850.6 


12:36:18.2 


+62 


15:51 


5.00 






7 


5 ± 0.9(7.5) 


170 + 7 


1.08 


0.001 


(107)MMJ123644+621147 


GN 850.28 


12:36:44.8 


+62 


11:47 


5.00 




^ 2.0 


1 


7 + 0.4(1.7) 


90 ± 24 


5.38 


0.031 


(108) MM 1123616+621517 


GN 850.4 


12:36:16.5 


+62 


15:17 


4.83 




^ 2.0 


5 


1 ± 1.0(4.9) 


58 + 7 


4.13 


0.031 


(109) MM J 123652+62 1353 


GN 850 30 


[2-36-52 7 


+62 


13:53 


4 75 




$ 2.0 


1 


8 ± 0.5(1.8) 


20 + 5 


1 89 


0.011 


/ 1 n> MM 1193737^69 1457 


nN 1 9nn i n 


j9.37.37 7 


+62 


14:57 


4 66 


2 


8 ± 0.6 (2.2) 




«: 6.0 








(33) MM J 123641 +621 130 


GN 1200 33 


12:36:41 6 


+62 


11:30 


4 6(1 


2 


6 ± 0.7(1.9) 




«: 3.0 








(7) MM J 123607 +620858 


GN 1200 7 


12-36-07 8 


+62 


08:58 


4 57 


3 


4 ± 0.7(2.8) 




^ 19.4 








(110) MM J 123608 +62 1433 


GN 850 23 


12-36-08 5 


+62 


14:33 


4 50 




^ 2.0 


7 


± 1.9(4.9) 


48 + 12 


1 4^ 


0.017 


(8) MM J 123702+62 1 950 


GN 1200.8 


12:37:02.6 


+62 


19:50 


4.42 


3 


3 ± 0.7(2.7) 




^ 10.7 








(13) MM J 123631 +621717 


GN 1200 13 


12-36-31 8 


+62 


17:17 


4 49 


3 


± 0.7 (2.2) 




^ 15.5 








(111) MM J 1 23622+62 1618 


GN 850 9 


12-36-22 2 


+62 


16:18 


4 33 




^ 2.0 


8 


9 ± 1.0(8.9) 


< 19.5 






(112) MM J 123709+6222 12 


GN 850 20 2 


12-37-09 9 


+62 


22:12 


4 99 




^ 2.7 


11 


.7 ± 2.2(9.9) 


160 ± 9 


6 45 


0.024 


(113) MM J123645+621447 


GN850.12 


12:36:45.1 


+62 


14:47 


4.16 




^ 2.0 


8 


6 ± 1.4(8.6) 


103 ± 6 


6.87 


0.038 


(21) MM J123720+621108 


GN 120021 


12:37:20.4 


+62 


11:08 


4. 14 


2 


8 ± 0.7(2.0) 




10.0 








(17) MM J 123735 +621 104 


GN 1200 17 


12-37-35 8 


+62 


11:04 


4 14 


2 


9 + 0.7 (2.1) 




«: 3.0 


42 + 7 


7 


061 


(16) MM 1 1 23723 +622050 


GN 1200.16 


12:37:23.8 


+62 


20:50 


4.00 


3 


3 ± 0.8 (2.3) 




^ 4.7 








(20)MMJ123749+621560 


GN 1200.20 


12:37:49.0 


+62 


16:00 


4.00 


2 


8 ± 0.7(2.1) 




< 23.9 








(114)MMJ123702+621304 


GN850.18 


12:37:02.4 


+62 


13:04 


4.00 




^ 2.0 


3.2 ± 0.6(3.2) 


21 + 6 


1.95 


0.013 


(115)MMJ123711+621245 




12:37:11.5 


+62 


12:45 


4.00 




^ 2.0 




^ 4.8 








(29)MMJ123711+621218 


GN 1200.29 


12:37:11.8 


+62 


12:18 


4.00 


2 


7 ± 0.6 (2.0) 




^ 4.0 








(116) MM 1123713+621200 


GN 850.21 


12:37:13.7 


+62 


12:00 


4.00 




^ 2.0 


5.7+ 1.2(4.9) 


34 + 6 


4.12 


0.040 


(27)MMJ123607+621138 


GN 1200.27/GN 850.2 


12:36:07.2 


+62 


11:38 


4.00 


2 


7 ± 0.7(2.0) 


16.2 ± 4.1 (12.1) 












3.5 ^ 


O" < 4 detections 


















(19) MM J 123645 +620824 


GN 1200.19 


12:36:44.9 


+62 


08:24 


3.85 


2 


8 ± 0.7(2.1) 




^ 7.3 








(42)MMJ123700+ 621455 


GN 1200.42/GN 850.32 


12:37:00.0 


+62:14:55 


3.83 


2 


1 ± 0.6 (1.4) 


5 


3 + 1.4 (3.8) 


53 ± 13 


5.04 


0.043 


(117)MMJ123618+621001 


GN850.5 


12:36:18.8 


+62:10:01 


3.83 




^ 1.9 


6 


7 ± 1.6 (5.2) 


50 ± 15 


4.14 


0.034 


(118) MM J '123638+621210 




12:36:38.4 


+62:12:10 


3.80 




2.0 




s; 5.9 








(119) MM J 123803+621446 




12:38:03.5 


+62:14:46 


3.75 




^ 2.1 




^ 2.5 








(23) MM J 1237 17+620801 


GN1200.23 


12:37:17.3 


+02 


08:01 


3.75 


3 


1 ± 0.8 (2.2) 




s; 3.0 


97 ± 19 


3.38 


0.015 


(28) MMJ12360S+6216S1 


GN1200.28 


12:36:05.6 


+62:16:51 


3.71 


2 


7 ± 0.7(2.0) 




16.9 








(32) MM] 1237 11 +621924 


GN 1200.32 


12:37:11.5 


+62:19:24 


3.62 


2 


6 ± 0.7(1.9) 




^ 10.9 








(120) MM J 123755+62 17 18 




12:37:55.3 


+62:17:18 


3.62 




^ 2.1 




^ 46.4 








(121) MM J 123656+621648 




12:36:56.8 


+62:16:48 


3.57 




^ 2.6 




^4.6 








(122) MM 1123621+621238 




12:36:21.2 


+62:12:38 


3.57 




^ 2.0 




«S U-5 








(123)MMJ123802+ 621503 




12:38:02.4 


+62:15:03 


3.50 




2.0 




15.1 








(14) MMJ123635+620715 


GN 1200.14 


12:36:35.5 


+62:07:15 


3.50 


3 


± 0.7(2.2) 




^ 19.5 


42 ± 13 


7.60 


0.073 



MAMBO/SCUBA associations in GOODS-N, and therefore 
are dealing with small number statistics, the sources are seen to 
span a similar range in 5 , 850 M m/S'i200nm and Ssso^m as existing 
samples. Most sources in all three samples exhibit a correlation 
between S , 850nm/S l i2oo M m and S^o^m- There are, however, a 
few outliers, although the uncertainties on some of these are 
substantial. 

G04 (contrary to E03) did not find a significant fraction of 
sources with low Sssonm/Si^ooum values, but there is no dis- 
crepancy between E03 and the G04+I05 sample compiled here. 



This is due to the deboosted Sssoum values and the wider and 
deeper coverage at 850 /im. As a result of the good agreement, 
we have combined all three samples in Fig. [8] in order to im- 
prove the statistics. The combined flux ratio distribution, shown 
as the empty histogram, is seen to be asymmetric around its 
peak at 5 , 850, lm /Si200nm ~ 2.5, with 12 out of 48 sources (25 
per cent) with S l 850 ( im/Si200Mm ^ 3 and a shoulder of low 
S , 850(im/<S'i200Mm values, with 18/48 (38 per cent) of the sources 
lying at 5850Mm/S'i2oo A1 m ^ 2. 
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Figure 8. Sgsuum versus the 5 8 50Mm/Sl200Mm ratio - with correspond- 
ing uncertainties, for the deboosted sources identified in both the SCUBA 
and MAMBO source catalogues of the GOODS-N field (red symbols), the 
sample of 16 ^3.5-cr SCUBA sources associated robustly with MAMBO 
in the Lockman Hole and ELAIS N2 fields (blue symbols) and the sam- 
ple of MAMBO sources with SCUBA photometry presented by E03 (green 
symbols). The open histogram shows the combined S , 850fim/-S'i200fim dis- 
tribution of all three samples, while the filled grey histogram shows the dis- 
tribution for the subset of radio-identified sources. 



5.2 Sources detected by MAMBO 

Next we expand our sample to include the MAMBO sources 
(20 sources - see Table [T) plus MAMBO sources which were de- 
tected at ^4ct in the combined map (GN 1200.21, .27, .29, .33). 
The total tally is 18 sources within the SCUBA region, with six 
sources outside (see bottom panel of Fig. |7J. In cases where there 
was no robust SCUBA association, we adopted the 3-a limit in the 
850-^tm map at the position of the MAMBO source. The result- 
ing plot of S850 M m versus Ssso^m/S^oo^m is shown in Fig. [5Ji. 
The uneven noise properties of the SCUBA map are reflected in 
the large range of flux density ratios exhibited by the MAMBO 
sample. A handful of MAMBO sources in the shallowest parts of 
the SCUBA map have poorly constrained flux ratios (^10). On 
the opposite side of the plot, we find three sources with tight up- 
per limits on fiasc^m, implying Ssso^/S^oo^m S 2. A num- 
ber of such sources are also found when repeating the analysis 
for the ELAIS N2 and Lockman Hole fields (S02+G04). For these 
two fields the SCUBA and MAMBO maps were well-matched in 
depth and, as a result, no sources with very high flux density ra- 
tios are found. Combining the S02+G04 results with those from 
GOODS-N, we obtain the open histogram shown in Fig. [9^. Ac- 
counting for the upper limits in the distribution, we find that 20- 
35 (37-65 per cent) of the 54 sources in the combined distribu- 
tion lie at S , 850Aim/5'i200Mm ^ 2 while 9-18 (17-33 per cent) have 

5 , 850^m/<S'l200M 
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Figure 9. a) Sssoum versus Sg^^m/S^oo^m for the entire ^ 4 " CT 
MAMBO sample and MAMBO sources detected at ^4tr in the combined 
map (red symbols). For sources with no robust SCUBA association we used 
the 3-<r flux density limit at the MAMBO position in the SCUBA map. 
The sample of MAMBO sources obtained from the Lockman Hole and 
ELAIS N2 fields (S02+G04) with SCUBA associations or upper limits on 
S'850fim are also shown (blue symbols). The flux ratio distribution obtained 
by combining the two samples is given by the open histogram, while the dis- 
tribution of radio-identified sources is shown as the grey histogram, b) The 
opposite experiment from a): SgSOum versus Sssoiim/ S^oo^im for ^ e 
sample of SCUBA sources in GOODS-N and SCUBA sources de- 

tected at ^4 a in the combined map (red symbols). SCUBA sources from 
S02+G04 with MAMBO counterparts or S^oOfim upper limits (3<r) are 
shown as blue symbols. The open histogram shows the flux ratio distribu- 
tion for the combined GOODS-N and S02+G04 samples while the grey his- 
togram shows the distribution for only radio identified sources (see t]5-4i . 
In both a) and b), the dotted line illustrates the 3-cr detection limit in the 
deepest part of the SCUBA map (~0.9 mly - Pope et al. 2005), while the 
dashed line shows the value of SssOpm corresponding to the 3-cr MAMBO 
detection limit (~2.1 mly) at a given flux density ratio. 
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5.3 Sources detected by SCUBA 

As an important check, we performed the mirror-image experiment, 
deriving S850Mm/£i200/mi for the SCUBA sample. We included all 
SCUBA sources (22 sources - Pope et al. 2006) and sources 
which were detected at a in the combined map (GN 850.39, .23, 
.28, .30 and .40), yielding a total of 27 sources. As expected, most 
of the sources lie along the dashed line in Fig. |9p, which shows 
the flux density ratio corresponding to the 3-cr detection limit of 
the MAMBO map (~2. 1 mjy) for a given 850-/im flux density. 
The flux ratio distribution obtained by combining the GOODS- 
N sample with that of S02+G04 has a single dominant peak at 
S850Mm/S'i200Mm ~ 2 — 3. Taking into account the lower limits 
in the combined distribution, we find that 10-37 (11-39 per cent) 
of the 94 sources in the distribution lie at S'850/jm/>S'i200/jm ^ 2, 
while 26-74 (28-79 per cent) lie at S850 Mm /<5i200/mi ^ 3. 

A Kolmogorov-Smirnov (KS) test of the flux density ratio dis- 
tributions of the combined MAMBO and SCUBA samples in Fig. [9] 
yields a likelihood of ~46 per cent that identical distributions could 
differ by as much as that observecQ, i.e. the distributions are indis- 
tinguishable statistically. Nonetheless, we do note that 37-65 per 
cent of the sources in Fig.|9^ have Ssso^m/Siwo^m ^ 2 and 17- 
33 per cent have Ssso/im/Si^oo/im ^ 3, while the corresponding 
percentages in Fig.|9p are 11-39 and 28-79 per cent, respectively, 
i.e. the fractions are almost reversed between the two distributions. 

Due to the uneven noise properties of the SCUBA map in par- 
ticular, the non-negligible errors on the (sub)mm flux densities, and 
the small sample sizes, we are unable from Fig.|9]to make conclu- 
sive statements about any difference between the Ssso/zm/Si^oo/im 
distributions of MAMBO- and SCUBA-selected samples. How- 
ever, we will reinvestigate this same question using two more fruit- 
ful approaches in §[6]and|7] But before doing so we will finish this 
section by examining whether radio identifications are helpful here. 



5.4 850-/1200-^m flux ratio distribution of radio-bright and 
radio-dim (sub)mm sources 

For the SCUBA samples in the Lockman Hole and ELAIS N2 re- 
gions we adopted the radio-identified SMGs in Ivison et al. (2002, 
2007). For the SCUBA sample in GOODS-N, we adopted the ra- 
dio identifications by Pope et al. (2005, 2006) and Biggs & Ivison 
(2006). To identify radio counterparts to our MAMBO sample in 
GOODS-N, we used the Biggs & Ivison (2006) radio catalogue to 
search for sources within 6arcsec of each MAMBO source. For 
every such association, the probability, P, of this happening by 
chance was calculated using the method of Downes et al. (1986). A 
radio identification was considered robust if P < 0.05. In TablesQ] 
and |21 respectively, we list radio identifications for the catalogues 
extracted from the MAMBO data alone, and from the combined 
SCUBA/MAMBO map. 

We find that only 9/30 (30 per cent) of the deboosted MAMBO 
sources have robust radio identifications, significantly lower than 
the typical identification fraction (65 per cent - e.g. Ivison et 
al. 2002) for bright SCUBA sources and radio data of similar 
depth (<ti.4ghz ~ 5 /iJy beam -1 ). Considering only the 1200- 
/im sources detected at J?4 a, the fraction increases although to no 
more than 35 per cent. If we instead consider the sample extracted 
at ^4 a from the combined map (Table O, we find 17/33 radio 



1 In performing the KS tests between the MAMBO and SCUBA samples 
we have assumed that the upper and lower limits are measured values. 



identifications (52 per cent), roughly consistent with that of bright 
SCUBA sources. 

From the combined samples in Fig.[8]we find that the fraction 
of MAMBO/SCUBA associations which also have radio identifi- 
cations is roughly constant (~60-100 per cent) across the range of 
flux ratios spanned by the combined sample. This is not surprising 
since we have only considered sources with high-S/N SCUBA and 
MAMBO detections. We therefore conclude that SMGs exhibit a 
broad range of Sssonm/Sisooum values, including quite low ratios, 
independent of whether or not they are identified at radio wave- 
lengths. 

Turning our attention to Fig. [9] and the flux ratio distributions 
of the MAMBO- and SCUBA-selected samples, and in particular 
the fraction identified at radio wavelengths as a function of flux 
density ratio (grey histograms), we find that at low flux density ra- 
tios (S85o,jm/Si2oonm S 2) the fraction of MAMBO sources with 
radio identifications is 0-60 per cent, while the corresponding frac- 
tion for the SCUBA-selected sample is 47-100 per cent. Given that 
this fraction is expected to drop at z £ 3, this is the kind of scenario 
we would expect to see if 1200-^im observations were selecting a 
significant population at those redshifts. However, these are small 
number statistics; a formal KS test does not support the notion that 
the two populations differ, giving a 24 per cent chance that the two 
distributions are drawn from the same parent distribution. 



6 STACKING ANALYSIS 

As we saw in Sj4] only five ^4-cr sources in the combined map 
of the GOODS-N region are robustly detected at both 850- and 
1200-/im. In this section we shall examine whether the remaining 
MAMBO (SCUBA) sources with no associations are detected at 
850 f-im (1200 /im) in a statistical sense. 

To this end, we have carried out a 1200-/im stacking anal- 
ysis of the 16 SCUBA sources in Table [2] with no MAMBO 
associations, but which were extracted from the combined 
SCUBA/MAMBO map at significance Js4 a. To avoid contamina- 
tion by bright MAMBO sources and their negative sidelobes we 
produce a CLEANed version of the MAMBO map in which all 
the ^3.5-0" deboosted sources (including their negative sidelobes) 
were masked, as well as a region 1 1 arcsec in radius around each 
one. The 1200-/im signal and noise values were extracted from 
the MAMBO map at the positions of the SCUBA sources, and the 
(variance-weighted) mean flux density and noise values were de- 
rived. This yielded a total of 15 SCUBA sources with a stacked 
flux density of (Si2oo M m} = 1-0 ± 0.1 mjy (S/N = 10). The aver- 
age Ssso/jm for these 15 SCUBA sources was 3.8 ± 0.2 mjy. As 
a check, we repeated the stacking analysis for the sample 
in Tableland found good agreement with the stacked flux density 
of the ^4-cr sample. As an extra check, we performed the stack- 
ing analysis with 10 per cent of the brightest and faintest of the 
1200-/iin flux density values removed from the stack, again finding 
agreement. 

In a similar manner, we derived the stacked 850-/^m flux den- 
sity of the 1 1 MAMBO sources in Table [2] with no SCUBA as- 
sociations, but which were recovered at ^4 a from the combined 
map. To this end, we used the CLEANed version of the SCUBA 
signal map (Pope et al. 2005) to ensure the result would not be af- 
fected by the negative sidelobes around each source. A 15-arcsec 
region around each ^3.5-<r source was also blanked out. The fi- 
nal stack contained all 11 sources and yielded a stacked flux of 
(S850nm) = 1.7 ± 0.6 mjy (S/N = 2.8). The average Ssso^m of 



© 0000 RAS, MNRAS 000, 000-000 



A MAMBO survey of the GOODS-N field 1 3 



the 11 MAMBO sources was 2.4 ± 0.2 mjy. We found that the 
stacked Sssoum f° r the sample was consistent with the ^4- 

a sample and that removing the 10 per cent brightest and faintest 
sources from the stack did not change the result. Our stacking anal- 
ysis is summarised in Table [3] where we have also listed the range 
in SssoAun/Siaoo/un implied by our analysis. 

To gauge the significance of these results, we performed 
Monte Carlo simulations in which the stacked signal was deter- 
mined in the same manner as above, except with randomised source 
positions. For the randomisation, we only considered positions 
within the MAMBO/SCUBA overlap region. Furthermore, since 
we wanted the simulations to represent similar noise properties to 
those in the original stacking analysis, each source was assigned 
a random position 20-80 arcsec from its original position. From 
10,000 such runs - each time with different positions - we found 
that the observed 1200-^tm signal occurs on ^0.001 per cent of oc- 
casions, while there is a 4 per cent chance of obtaining the observed 
850-^tm stacked signal. This confirms the high significance of the 
stacked 1200-^im signal and suggests the 850-^im signal could oc- 
cur by chance. 

It is important to realise that these stacked fluxes are likely 
to be biased low, since the MAMBO and SCUBA maps are un- 
likely to be perfectly aligned and since the positional uncertainties 
of (sub)mm sources at low significance levels are ~4 arcsec. As 
a result, we repeated the entire stacking analysis, this time using 
the peak flux within lOarcsec of each source in the stack, rather 
than the flux density at the exact position (see § [4] for the reason 
for choosing a 10-arcsec radius). This yielded a stacked flux of 
(Si2oo M m> = 1.8 ± 0.1 mjy (S/N = 18) for the 1200-Atm-blank 
SCUBA sources; for the 850-/im-blank MAMBO sources we found 
(SWm) = 4.0 ± 0.7 mjy (S/N = 5.7). 

While this method yields an apparently significant 850-/im de- 
tection of the 850-/im-blank MAMBO sources, selecting the peak 
within lOarcsec of each source will bias the stacked 850-^tm flux 
density to high values. We consider the estimated signal from this 
method to represent an upper limit, so the average 850-/im flux den- 
sity lies in the range (Sssofwn) — 1.7 — 4.0 mJy. The SCUBA sam- 
ple is detected significantly (^ 10 <r) at 1200 fim with both methods 
with (Si200Mm) — 1.0 — 1.8 mjy. 

Stacked fluxes obtained by selecting peaks within a cer- 
tain radius are prone to flux boosting. We therefore repeated the 
above analysis, this time deboosting the peaks using the Bayesian 
scheme described in § 13.2.31 This yielded a stacked flux of 
(SiaooMm) = 1.2 ± 0.1 mjy (S/N = 12) for the 1200-^m-blank 
SCUBA sources, while for the 850-/im-blank MAMBO sources we 
found (SssOMm) = 3.0 ± 0.7 mjy (S/N = 4.3). The net effect of de- 
boosting is to lower the stacked flux estimates, narrowing the range 
in (S , 850Mm)/<Si2oo M m} for both the MAMBO and SCUBA sam- 
ples (see Table[3j. 

In addition to the stacking analysis, we compared the distri- 
bution of 850- ^im S/N values at the positions of the 850-pim-blank 
MAMBO sources with that of the overall S/N distribution of the 
(CLEANed and blanked) SCUBA map, and similarly compared the 
1200-^m S/N distribution of the 1200-/im-blank SCUBA sources 
with the overall S/N distribution of the (CLEANed and blanked) 
MAMBO map (only considering the region overlapping with the 
SCUBA map). Comparing the full S/N distributions is likely to 
be a more robust method than calculating the stacked flux densi- 
ties, which only uses the first moment. We performed this anal- 
ysis for the 1200-/im-blank SCUBA sources and 850-^tm-blank 
MAMBO sources which had been extracted at ^2.5 a from the 
SCUBA and MAMBO maps, respectively. We stress that by ' 1200- 



fim- (850-/im-)blank\ we mean sources which did not have a 
>2.5-<r MAMBO (SCUBA) association within 10 arcsec of its po- 
sition in the combined map. The resulting distributions are shown 
in Fig. [10] The 1200-yum S/N distribution of the 1200-/im-blank 
SCUBA sources is clearly skewed towards positive values rela- 
tive to the map. A KS test gives a 0.013 per cent likelihood that 
the two distributions are identical. For the 850-/im-blank MAMBO 
sources, the 850-^im S/N distribution is not strongly inconsistent 
with that of the overall SCUBA map (KS test probability, 4.5 per 
cent). 

The statistical detection of the 1200-/im-blank SCUBA 
sources is consistent with these sources lying just below the detec- 
tion threshold at 1200-/im, but being otherwise identical to the de- 
tected sources. The non-detection of the 850-/im-blank MAMBO 
sources - except when non-deboosted peak values are used in the 
stack - suggests that these sources lie well below the detection 
threshold, being much fainter than typical SMGs. These results are 
also reflected in the average Ssso/Mn/Siaoo/Mn values. For the lat- 
ter, £ , 850Mm/S'i200Mm £ 1.7, while the 1200-^m-blank SCUBA 
sources have Sssoitm/Suioopm values consistent with SMGs which 
have robust detections at both 850 and 1200 /im. Our findings 
therefore support the notion that a significant fraction of MAMBO 
sources selected from blank-field surveys are SDOs, with cooler 
1200-850 fim colours. 



7 SIMULATIONS OF THE ASSOCIATION FRACTIONS 

Is the observed fraction of SCUBA sources with MAMBO asso- 
ciations (and vice versa) consistent with what we would expect 
given the properties of the maps? Since the noise varies signifi- 
cantly across the SCUBA map, we must address this question by 
means of careful Monte Carlo simulations of the associations, tak- 
ing into account the varying depths of the SCUBA and MAMBO 
maps. 

To this end, we took the ^2.5-er SCUBA and MAMBO cata- 
logues and searched for counterparts within 10 arcsec of each other. 
In total, 28 such SCUBA/MAMBO associations were found in 
this manner: 63 per cent of the SCUBA sources had 1200- 

/im associations; 30 per cent of the ^4-<r MAMBO sources had 
850-/^m associations. From the 28 SCUBA/MAMBO associations, 
we constructed a <S , 850Mm/S , i200/jm distribution (see Fig.fTTb), tak- 
ing into account the errors (assuming Gaussianity). If SCUBA and 
MAMBO sources - even at faint flux density levels - trace the same 
galaxy population, we would expect the SssOfim/Smo/jm distribu- 
tion for the ^2.5-ct samples to be representative for all SCUBA and 
MAMBO sources at all significance levels. This is the null hypoth- 
esis we want to test. 

First, we examined whether the observed fraction of SCUBA 
sources with MAMBO associations is consistent with what we 
would expect. We ran a Monte Carlo simulation in which we as- 
signed a random flux density - drawn from the resampled 850-/im 
flux density distributions of ^4-a SCUBA sources (Fig. II lb ) - 
to each member of the sample of ^4-cr SCUBA sources detected 
with MAMBO at >2.5 a. Then a Sssofim/Siwopm value was ran- 
domly selected from the SssoMm/S^oOfim distribution, and the cor- 
responding S'i200fim was calculated. 

Next, we identified the regions of the SCUBA map where each 
850-^im flux density value could be detected at ^4<r, given the 
observed noise. Since the noise varies across the SCUBA map, so 
the suitable regions vary for each flux value. We then drew random 
noise values, <7i200fim , from the noise distribution of the parts of the 
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Figure 10. Left: Distribution of 12()0-/im S/N values for the SCUBA sources (hashed histogram and left j/-axis) compared to the 1200-/im S/N 

distribution of the CLEANed and blanked MAMBO map (open histogram and right y-axis). The former is skewed towards positive values and a KS test yields 
a 0.013 per cent likelihood that the two distributions are identical. Also shown for comparison is the distribution of the 1200-/im S/N values of the ^2.5-cr 
MAMBO sources (grey hashed histogram), normalised to the histogram for the ^2.5-cr SCUBA sources. Right: Distribution of the 850-/^m S/N values of 
the JS2.5-CT MAMBO sources (hashed histogram and left y-axis), compared to the 850-/im S/N distribution of the CLEANed and blanked SCUBA map (open 
histogram and right y-axis). The two distributions are drawn from the same parent distribution with 4.5 per cent likelihood. For comparison, the 1200-/^m S/N 
distribution of the Jj2.5-<r SCUBA sources is also shown (grey hashed histogram). 



Table 3. Summary of stacking results for the 1200-/xm-blank SCUBA sources and the 850-/im-blank MAMBO sources. Both the signals stacked using the 
flux densities at the (sub)mm centroids (biased low) and the peak flux densities within a 10-arcsec aperture (biased high) are listed. The last column lists the 
stacked signals obtained when deboosting the peak fluxes. The percentages in parentheses are the likelihoods of obtaining the stacked fluxes by chance. 





1200 
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1.0 ± 0.1 (< 10 -3 ) 
3.8 ± 0.2 
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MAMBO map which overlapped with the suitable SCUBA regions. 
In doing this, we have accounted for the varying noise across the 
maps. An example of this is shown in Fig.ll lfc. where we show the 
regions of the MAMBO noise map which overlap with the regions 
in the SCUBA map where a flux density of 7. 1 mJy is detectable 
at ^4<r. The noise distribution of this region is shown as the grey 
histogram in 

Fig.HB- 

Finally, the fraction of associations was calculated, based on 
how many sources had S^oo^m/a^oo^m ^ 2.5. We ran this ex- 
periment 10,000 times, building up a distribution of association 
fractions, where by 'association fraction' we mean the fraction of 
SCUBA sources with associations in the MAMBO image. From 
the resulting distribution, shown in Fig. II lfe . we see that that the 
observed 1200-/im association fraction (63 per cent) is fully consis- 
tent (within 1 a) with the simulated distribution of association frac- 
tions (median, 75 per cent). We ran two additional simulations, sim- 



ilar to that described above, but using the ^3.5- and SCUBA 
samples. In both cases, we found that the observed association frac- 
tion is consistent with the simulations. The observed fraction of 
SCUBA sources recovered at 1200 /j,m is expected, then, given the 
noise properties of the maps. 

Can the same be said for the 850-/iin recovery fraction of 
MAMBO sources? To test this, we carried out the mirror-image 
experiment as above, i.e. random 1200-^tm fluxes were assigned to 
the sample of ^4-<r MAMBO sources with ^2.5-a SCUBA detec- 
tions, using the resampled 1200-/im flux distribution (Fig. 
[TTh . Using the same 850-/1200-/im distribution as above (i.e. de- 
rived from the ^2.5-cr detections) the 850-/im flux densities were 
determined for the sample. In a similar manner as before, the re- 
gions of the MAMBO maps in which each 1200-/im flux density 
could be detected at ^4<r were identified and 850-/im noise val- 
ues were randomly extracted from the the noise distributions of the 
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MAMBO r.m.s. noise [mJy/beam] ID fraction S,„„ [mJy] 




111 <arcs = c > SCUBA r.m.s. noise [mJy/beam] ID fraction 

Figure 11. Distributions from our Monte Carlo simulations, described in §7] a) SsSOum/SwoOiMn distribution derived by resampling the flux density ratio 
distribution of the 28 SCUBA/MAMBO associations found as described in § [7] taking into account the flux density errors, b) The 850-/xm flux density 
distribution obtained by resampling the flux distribution of the ^A-o SCUBA sources. The vertical dotted line represents a flux density of 7. 1 mJy drawn from 
this distribution, c) The grey-scale regions show the parts of the MAMBO noise map which correspond to regions in the SCUBA map where a flux density 
of 7. 1 mJy can be detected at cr. This serves as an example of how we have taken the inhomogeneous noise distribution of the maps into account, d) The 
open histogram shows the noise distribution of the entire MAMBO map, while the noise distribution of the regions selected in c) are shown as a filled, grey 
histogram, e) The resulting 1200-/im association fraction distribution (normalised to unity) for the SCUBA sample. The vertical solid line indicates the 

median of the distribution (75 per cent), while the dotted line shows the actual observed fraction (63 per cent). The solid curve represents a single Gaussian 
fitted to the distribution, f) The flux density distribution of the Jj4-cr MAMBO sample. The vertical dotted line represents a randomly chosen flux value of 
2.8 mJy. g) The grey-scale regions show parts of the SCUBA noise map which correspond to regions in the MAMBO map where a flux density of 2.8 mJy 
can be detected at cr. h) The noise distribution of the entire SCUBA map (open histogram), and the noise distribution corresponding to the noise regions 
shown in g) (filled, grey histogram), i) The 850-/^m association fraction distribution (normalised to unity) for the Js4-<r MAMBO sample. The vertical solid 
line indicates the median of the distribution (74 per cent) and the dotted line shows the actual observed fraction (30 per cent). 
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matching regions in the SCUBA noise map. Finally, the fractions 
of associations were calculated. The resulting distribution of 850- 
/im association fractions for the MAMBO sources is shown 
in Fig.ll 11 and it is seen that the observed 850-pim association frac- 
tion (30 per cent) of MAMBO sources falls completely outside the 
simulated distribution of association fractions, deviating from the 
mean (74 per cent) by 4 a. 

The conclusion of this analysis is that, provided we are pre- 
pared to use the sample of lower significance peaks, we can re- 
ject the null hypothesis for MAMBO sources. In other words there 
is evidence for a population of 1200-^im sources which does not 
share the distribution of flux density ratios with sources found 
jointly at 850 and 1200 (jm. The non-trivial distribution shown in 
Figs. Illfe and i underscores the importance of performing these 
careful Monte Carlo tests in order to be confident of this result. 
In order to explain the low fraction of associations, given the sig- 
nal and noise properties of the two maps, the 850-/1200-/im flux 
density distribution for the MAMBO population would have to be 
skewed towards lower values. This argues in favour of a significant 
fraction of them being SDOs, at higher redshifts and/or with cooler 
dust. 



8 DISCUSSION 

The spectroscopic redshift survey of bright, radio-identified SMGs 
by Chapman et al. (2003, 2005) located sources out to z ~ 3.6, 
with an interquartile range, 1.7 z ^ 2.8. 105 argued that of the 
bright SMG population (Ssso/im 5 mjy), probably no more than 
~10 per cent could be at z > 3.5, given that ~80 per cent have ra- 
dio counterparts and ~10 per cent were likely to be spurious. Pope 
et al. (2006) reached a similar conclusion, estimating that <14 per 
cent of SMGs reside at z > 4. Modeling the radio/mm/far-infrared 
colours of 120 SMGs from SHADES and adopting priors for the 
redshift probability of their radio-undetected sources, Aretxaga et 
al. (2007) argued that more than half of the bright SMG population 
lies in the range 1.6 & z ^ 3.4, with little room for a z > 4 pop- 
ulation. Thus, the evidence for a significant population of SMGs at 
very high redshifts appears slim, but has not been ruled out com- 
pletely. 

If we consider sources strictly selected at 850 fim, seven have 
spectroscopically confirmed redshifts in the range 3 < z < 3.6 
(Ledlow et al. 2002; Chapman et al. 2003, 2005). A strong candi- 
date for a z ^> 3 SMG was reported by Knudsen, Kneib & Egami 
(2006) who inferred a likely spectroscopic redshift of ~ 4 for 
SMM J16359+66130. Apart from this source, however, only one 
other convincing candidate for a z > 4 SMG has been presented 
in the literature at the time of writing, namely GN 850.10 (Wang 
et al. 2004; Pope et al. 2005). This source was deemed to reside at 
z ~ 4 — 6, based on its near- and mid-infrared colours, which were 
facilitated by its accurate location at 890-/im with the SMA (Wang 
et al. 2007), at 1.1 mm with the IRAM PdBI and in the radio with 
the VLA (Biggs & Ivison 2006; Dannerbauer et al. 2008). Younger 
et al. (2008) presented SMA observations of LH 850.2 - one of 
the brightest 850-/im (and 1200-/^m) sources in the Lockman Hole 
(G04; Coppin et al. 2006) - and derived an optical/infrared photo- 
metric redshift of z ~ 3.3, at the high-end of (but within) the Chap- 
man et al. redshift distribution. Aretxaga et al. (2007) presented a 
handful of SMGs from SHADES with photometric redshifts z > 4, 
but with very large individual uncertainties. 

Turning our attention to sources strictly selected at mm wave- 
lengths, we find a larger number of possible z > 4 candidates. Us- 




Z 



Figure 12. Green, blue and red filled circles represent S , 850fim/'5'i200Mm 
for sources identified robustly at 850 and 1200 fim by E03, S02+G04 and 
in this work, respectively. The filled triangles are the seven sources from 
Younger et al. (2007), where we have converted the measured 890- and 
ll()0-/tm flux densities to 850- and 1200-/im flux densities assuming an 
optically thin grey-body with /3 = 1.5. Sources which have robust radio 
counterparts have been circled. Where available, we have placed sources at 
their spectroscopic (Chapman et al. 2005) or photometric redshifts (Aretx- 
aga et al. 2007; Pope et al. 2006; Wall, Pope & Scott 2008), otherwise 
the sources have been placed at fixed redshifts below z = 0.3 for clar- 
ity. Sources with spectroscopic redshifts are shown as large symbols while 
sources with only a photometric redshift - or lacking a redshift estimate 
entirely - are shown as small symbols. The grey shaded region corresponds 
to the range of Ss5o M m/5i200Mm spanned by grey-bodies with f) = 1.0 
and T d = 10 K (lower limit) and f} = 2.0 and T d = 70 K (upper limit). 
Large red squares and stars represent the range in SgsOtim/ Si200/j.m de- 
rived from the (centroid and peak) stacking analysis for the 850-^m-blank 
MAMBO sample and the 1200-/xm-blank SCUBA sample, respectively. 
The open upward- and downward-pointing triangles represent the lower and 
upper limits on the flux density ratio, respectively, obtained when using de- 
boosted peak values in the stacking analysis (see §[6). 



ing IRAM PdBI to accurately locate their positions, Dannerbauer 
et al. (2002, 2004) presented three examples of bright (Si200/«n Z 
3.5 mjy) MAMBO sources which, based on their radio, (sub)mm 
and near-infrared properties, were argued to lie at z > 4. The non- 
detection of these sources with Spitzer/IKS spectroscopy (Valiante 
et al. 2007) is consistent with this claim. Younger et al. (2007) pre- 
sented 890-/im interferometric SMA observations of seven very 
bright (Siiooum ^ 7.6 mjy) AzTEC sources selected at llOO/im 
and argued that the five radio-dim (Siagh* ^ 41 /ijy) sources 
in their sample were likely to reside at higher redshifts than the 
radio-identified sources due to their systematically higher submm- 
to-radio flux ratios, lower IRAC 3.6-8.0-/im flux densities and non- 
detections at 24 jj,m. 

Common to all of these z > 4 candidates, whether they are 
SCUBA-, MAMBO- or AzTEC -selected sources, is that they are 
very bright (i.e. luminous) and have been detected at both submm 
and mm wavelengths. The trend that the high-redshift candidates 
tend to be very bright was first noted by Ivison et al. (2002), who 
from deep radio imaging of a large sample of SMGs found evi- 
dence that very bright (Sssopm Z 8 mjy) SCUBA sources tend to 
have higher submm-to-radio flux density ratios than less luminous 
SMGs, indicating that they are at higher redshifts (Carilli & Yun 
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1999) and that strong luminosity evolution may have taken place 
(see also Pope et al. 2006). 

From the above summary (and following § [7j, one wonders 
whether submm- and mm-selected sources are drawn from iden- 
tical populations? Does the preponderance of z > 4 candidates 
amongst the mm-selected sources reflect a real difference between 
the two populations, or is this merely a coincidence? The multi- 
wavelength properties of mm-selected sources have been studied 
far less thoroughly than SMGs, because they were discovered more 
recently and are often thought to be mere mm-wave analogues to 
SMGs. In particular, no systematic spectroscopic surveys have been 
undertaken of mm-selected sources and their redshift distribution is 
therefore not known. 

A related question is: why are all of the z > 4 candidates 
amongst the brightest (sub)mm sources known? Is this merely a 
selection effect, reflecting the fact that we can currently detect only 
the brightest sources with (sub)mm and radio interferometers, or 
does it reflect a genuine evolutionary effect (Wall, Pope & Scott 
2007)? 

In § |6] we showed that while the 1200-/im-blank SCUBA 
sources are statistically detected at 1200 /im (S/N ~ 10), the 850- 
/im-blank MAMBO sources are not robustly detected at 850 /im 
(S/N ~ 2.8). In Fig.[72lwe plot Ssso^m/ 81200^ for these two sam- 
ples, using the stacking results listed in Table[3] The ranges of pos- 
sible values for the average Ssso^m/S^ooMm ratio for the two sam- 
ples only barely overlaps, even when accounting for the high and 
low biases introduced when using the centroid and peak values, re- 
spectively (see Sj6]l. If we deboost the peak fluxes which go into the 
stack, the allowed ranges do not overlap at all. On average, there- 
fore, the 850-/im-blank MAMBO sample has significantly lower 
S , 850/jm/S'i200Mm ratios than the 1200-/im-blank SCUBA sample. 

The strongest support for this result was given in § UJ based 
on Monte Carlo simulations of the association fractions in the two 
maps. Since we are working in the low-S/N regime, the fraction of 
sources recovered in two different maps will never be unity, so it 
is crucial to simulate the experiments. Careful simulations are re- 
quired because of the inhomogeneous noise in the maps. The main 
outcome of our Monte Carlo simulations is that while the fraction 
of SCUBA sources detected with MAMBO is consistent with the 
S , 850(jm/<S'i200Mm distribution of sources found in both maps, there 
are more MAMBO sources not detected by SCUBA than would be 
expected. Since the simulations account for the noise properties of 
the two maps, we were able to estimate the fraction of MAMBO 
sources which are undetected by SCUBA (compared with the ex- 
pectation if the populations were the same, but in the presence of 
the observed noise properties) and found this to be ~40 per cent. 

Thus, we have presented strong evidence which suggests that 
- while there is substantial overlap between the 850-/im- and 1200- 
/im- selected sources in GOODS-N - the two populations are not 
identical. A significant fraction of the 1200-/im-selected sources 
are unaccounted for by the 850-/^m data. Due to the statistical na- 
ture of our analysis and the careful way in which we have taken into 
account the uneven depths and noise properties of the two maps of 
the GOODS-N region, we are confident that this is a general re- 
sult that reflects the existence of two distinct, albeit overlapping, 
populations. 

We still need to establish whether SDOs are genuinely at z > 4 
or whether they have cooler dust temperatures at redshifts typical 
of radio-identified SMGs. Due to the degeneracy between dust tem- 
perature and redshift (Blain 1999), and without direct spectroscopic 
redshifts, it is not trivial to determine which of the two explanations 
is appropriate. 



We note from Fig. [T2] that if the typical SEDs of SDOs are 
given by T c \ — 10 K and /3 = 1.0, which is believed to be ex- 
treme not only compared to starbursts at the present day but also 
to SMGs (Kovacs et al. 2006), then their most likely average flux 
ratio (S , 850Mm/Si200Mm ~ 1.2) would indicate an average red- 
shift of z ~ 1.3. One such source is seen in Fig. [12] namely 
Lock850.27/LH 1200.7 (G04; 105) which has a flux density ratio 
of ~ 1 and lies at z = 1.21. Given the general lack of spectro- 
scopic redshifts for SMG samples, further examples of such sources 
may exist, making them viable candidates for SDOs. Adopting 
Td = 20 K and /3 = 1.5, close to the coolest dust temperature mea- 
sured for SMGs so far (Kovacs et al. 2006), results in an average 
redshift of z ~ 4.6. If we assume that SDOs have SED proper- 
ties similar to those of radio-identified SMGs, i.e. Td = 35 K and 
f3 = 1.5 (Kovacs et al. 2006), we find their average redshift to be 
~ 9, in the regime where the fc-correction starts to become positive 
(Blain & Longair 1993). 

In the first of the three scenarios outlined above, we note that 
the angular sizes of cool (Td S 25 K) SMGs at z < 2 should scale 
roughly as 9 ~ (Sssonm/mJy) 1 / 2 (Kaviani et al. 2003) which 
would imply that SDOs have typical sizes of 9 ~ 1.3arcsec (for 
S850nm = 1.7 mjy - see Table[3}, corresponding to a physical di- 
ameter of 11 kpc at z — 1.5. Thus if SDOs reside at z £ 2 they 
would have larger physical sizes than those observed for SMGs 
(~4kpc - Tacconi et al. 2006; Biggs & Ivison 2008), but similar 
angular sizes (~1 arcsec). If SDOs predominantly lie at the same 
redshifts as SMGs (z ~ 2.5) but have cooler dust temperatures 
(Td S 25 K), their angular sizes would be larger by a factor of 
two (Kaviani et al. 2003). Thus, if SDOs are in fact very cool sys- 
tems with a redshift distribution similar to SMGs (or peaking at 
lower redshift), we would expect them to be quite extended in high- 
resolution radio and mm images. 

If some SDOs, on the other hand, lie in the range z ~ 4 — 10 
and have typical dust temperatures in the range Td ~ 20 — 40 K, the 
expected angular and physical sizes are 9 ~ 1.5 — 2.3 arcsec and 
6.4-18.0 kpc, respectively, inconsistent with the typical sizes mea- 
sured for SMGs. Thus, if the Kaviani et al. (2003) results are appli- 
cable, it would suggest that SDOs are very extended systems over 
a broad range of plausible dust temperatures and/or redshifts. We 
note that high-resolution 1.3-mm interferometry of SDOs should 
be able to test this and may provide a useful way of discriminat- 
ing between typical SMGs - generally compact at mm wavelengths 
(^1 arcsec - Tacconi et al. 2006; Younger et al. 2007) and SDOs, 
which our analysis here suggest may be larger. 

Observational evidence for such extended systems was re- 
cently uncovered by Daddi et al. (2008), who demonstrated the 
existence of gas-rich disks at z = 1.5 with physical sizes 2-3 
times those of SMGs, as implied by CO observations as well as 
rest-frame UV light. Some of these BzK galaxies - so called for 
the way they are selected (Daddi et al. 2004) - have been detected 
at very faint flux levels at \2QQ-pm. (~1.5 mjy - Dannerbauer et al. 
2005) and have similar star-formation efficiencies to those of local 
spirals, i.e. an order of magnitude smaller than SMGs. They could 
be candidates for the cold, extended 'cirrus' dust models proposed 
by Efstathiou & Rowan-Robinson (2003). Certainly, these prop- 
erties seem to fit with those of SDOs, opening up the interesting 
possibility of an overlap between the two populations, with SDOs 
making up some fraction of the near-infrared-selected population 
at 2 ~ 1.5 - 2.4. 

However, it seems clear from a growing number of in-depth 
studies, that some of the brightest mm- (and submm-)selected 
sources are almost certainly at z > 4 (Dannerbauer et al. 2002; 
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Wang et al. 2004; Younger et al. 2007). The same may be true 
for many of the SDOs, although we stress that one cannot rule 
out a scenario in which SDOs are predominantly cool, low-redshift 
sources. 

Adopting the high-redshift scenario (z ~ 5) for SDOs, and 
assuming T d = 20 K, /3 = 1.5 and Sfesojun = 1-7 — 4.0 mjy 
(Table [3}, we find far-infrared luminosities and dust masses of 
Lfir = (1-2) x 10 12 Lq and M d = (3 - 8) x 10 9 M Q , 
respectively. In this scenario, SDOs are luminous and massive sys- 
tems, despite being relatively fainter at submm wavelengths. The 
presence of such systems at z ~ 5 (corresponding to ~1 Gyr after 
the Big Bang) would have implications for our understanding of 
galaxy formation and evolution. Due to their large masses and high 
redshifts, SDOs and SMGs at z > 4 would be potential candidates 
for galaxies caught in the act of collapsing early in the Universe's 
history (Eggen, Lynden-Bell & Sandage 1962). 



9 SUMMARY 

We have undertaken the first deep and uniform 1200-/im survey of 
the GOODS-N region, identifying 30 sources (after flux deboost- 
ing) in our image. By combining our 1200-/im MAMBO map with 
the existing 850-/im SCUBA data, we have extracted a robust sam- 
ple of 33 sources (sub)mm sources detected at a combined S/N 4. 

The principal question we wanted to address in this paper was 
whether the 850- and 1200-^im selected sources constitute iden- 
tical, or merely overlapping, populations. We performed three in- 
dependent statistical analyses of the 850- and 1200-/im maps and 
their corresponding source catalogues. 

The first - a simple comparison of the 850-/1 200- flux den- 
sity ratio distributions for the 850-/im- and 1200-/im-selected sam- 
ples, suggests that a larger fraction of MAMBO sources have low 
values of <S , 850fim/S'i200f 1 m compared to the SCUBA sources. This 
method suffers, however, from small sample sizes and the effects of 
maps with uneven noise properties, and so we were unable to see a 
significant difference between the two populations. 

However, using the other two methods - both of which take 
into account the uneven noise properties - we do find evidence 
to suggest that the source populations selected by MAMBO and 
SCUBA are not identical. The strongest evidence comes from our 
Monte Carlo analysis of the fractions of SCUBA galaxies with 
MAMBO associations, and vice versa, which showed that while the 
fraction of SCUBA sources with MAMBO associations was con- 
sistent with the map properties and the SssOfim/ Si200nm distribu- 
tion of sources robustly detected at both wavelengths, the fraction 
of MAMBO sources with SCUBA associations was significantly 
lower than expected. In fact, we found that about 40 per cent of 
the MAMBO sources were not recovered at 850 fim, thus lending 
strong evidence to the notion that the two populations are not iden- 
tical. 

We have argued that the average Sssopm/Siaooian °f SDOs is 
significantly lower than the bulk of the SMG population, suggesting 
that these sources are either very cool or lie at higher redshifts. It 
would require extremely cold and unusual (Td ~ 10 K and (3 ~ 
1) far-infrared/mm SEDs to explain the presence of SDOs at low 
redshifts (z < 2). Nonetheless, examples of such sources do exist. 
We hint at a link between such cool, extended mm-selected sources 
and the near-infrared-selected population of BzK galaxies, which 
are known to be gas-rich yet have star-formation rates an order of 
magnitude lower than those of SMGs. 

If we adopt more realistic SED parameters (Td ~ 20 K and 



~ 1.5) for SDOs, it would imply an average redshift of z 4, 
i.e. beyond the range probed by spectroscopic surveys so far. If this 
is the case, it is tempting to speculate that SDOs could be potential 
candidates for galaxies collapsing very early on in the Universe's 
history. 

It should be stressed, however, that our analysis merely shows 
that SDOs are at at generally higher redshifts (or are cooler) than 
SMGs. Extensive comparison with detailed models will be required 
in order to determine what range of redshift distributions are con- 
sistent with the available data. This was beyond the scope of the 
present study. This task will also be made much easier with addi- 
tional information on these sources, particularly spectroscopic or 
optical/IR photometric redshifts. However, this first requires find- 
ing an identification of these SDOs at other wavelengths. 

Deep interferometric observations at (sub)mm wavelengths 
sensitive enough to probe the faint flux levels of SDOs, are likely 
to become one of the most promising ways of furthering our un- 
derstanding of SDOs. Such observations would allow us to iden- 
tify near-/mid-infrared and optical host galaxies with which to 
constrain their redshifts. Due to the limited sensitivity of current 
(sub)mm interferometers, however, such studies have been limited 
to the brightest SMGs. With the new generation of wide-band re- 
ceivers, and with ALMA on the horizon, studies of larger sam- 
ples of SDOs are becoming feasible. Furthermore, the wide, deep, 
multi-colour surveys planned with SCUBA-2 and Herschel (and 
on a somewhat longer time-scale, with CCAT) may provide use- 
ful far-IR/submm photometric redshifts for huge samples of SDOs 
and SMGs (Hughes et al. 2002; Aretxaga et al. 2003). Armed with 
these redshift constraints, heterodyne receivers with extremely high 
bandwidth could be employed to search for redshifted CO and/or 
C II line emission. For a complete census of obscured star formation 
and AGN activity across cosmic time, we must determine the range 
of SED properties for (sub)mm-selected sources, along with their 
true redshift distribution. Our evidence that mm-selected galaxies 
may extend beyond z = 4 suggests that this may be an efficient 
way to select and study massive galaxy formation at the earliest 
times. 
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